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Estimated value of deck shear-lag factor ''r" is included in the theory. 
The longitudinal stress distribution is obtained using the developed 
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reasonable. The difference between the results at the end of deckhouse 
is more distinguishable. 
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INTRODUCTION 


In 1953, H. H. Bleich published a paper in the "Journal of Applied 
Mechanics" entitled "Non-Linear Distribution of Bending Stresses Due to 
Distortion of the Cross Section". In this paper he derived a viable 
analytical solution to the problem of hull-deckhouse interaction. 

Basically, he considered the hull and the deckhouse as separate 
beams which are forced to act together by shearing forces and by verti- 
cal forces resisting relative displacements of the two beams. The case 
of constant cross section of the beam is treated, and it is assumed 
that cri hypothesis is valid for the deckhouse and hull separately. 
For two types of loading he considered (constant moment loading and 
equally distributed load), solutions were in a qualitative agreement 
with the test results at midship. As one moves away from amidships or 
the center of the deckhouse structure, solutions were departing from 
the reality. 

It is proposed that including shear effects into Bleich's original 
two-beam deckhouse theory, it is possible to improve it, particularly at 
deckhouse ends. In order to confirm or disaprove this hypothesis, new 
theory which includes shear effects is developed, and the results are 
compared at midship and at the end of dette. 

The theorem of Minimum Potential Energy is used in the analysis. 
The variational procedure established two coupled six order differential 
equation systems and the natural boundary conditions. The boundary value 


problem is solved for the displacement and then the stress distributions. 
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FUNDAMENTAL ASSUMPTIONS 


The two components of the combination of hull and deckhouse are 
treated as if each were a Navier beam which are forced to act together 
by shearing forces and by vertical forces. 

It is assumed that the stiffness of bulkheads or deck beams 
resisting relative vertical displacement of: the deckhouse is constant 
for the full length of the deckhouse, the magnitude of the stiffness 
being given by a spring constant K. In reality, the deck stiffness K 
vary along the deckhouse length due to the presence of structural 
bulkheads. 

The possibility of having different materials for hull and deck- 
house is not considered in this analysis. i 


Shear deformation is accounted for in the sides of both hull 


and deckhouse. For this purpose, the following equations are employed. 


Py 
a tt + 
iin ea nee | (1) 
lw 
Ez (2) 
== - 11 RE 
M, = -E LO, tE G? 
2w 


The longitudinal stress in the deck at the junction and the 
longitudinal stress in the deck-edge may differ. Bleich's theory 
neglects this. To include the shear-lag effect in the deck-edge and 
deckhouse side, the following assumption is used, 


0. =Yr0U (3) 


where r is the ratio of longitudinal deck stress at the junction to 


the stress at deck-edge. 
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Precise determination of "r" would require two-dimensional 


elastic analysis of the response of all plate elements of the section. 
For the present purpose, it is assumed that "r" has the same value it 
would have without the deckhouse, determined by a box-girder analysis 
of the hull alone, such as in Reference (2). According to this analysis, 


for a sinusoidal bending moment: 


r = u tnt 2 an a La T 75 (4) 
2 mb £ 2 4 L L L 
cosh -7 


It is pointed out in Reference (4) by Shade, for a bending 
moment which is constant over the length of the deckhouse "r" should 


be taken as unity. 
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CHAPTER I - ANALYSIS OF A SIMPLIFIED TWO CELL STRUCTURE 


Consider the problem of two separate beams forced to act together 
by horizontal shear forces and vertical forces acting at the junction of 
hull and deckhouse, (Figures 1 and 2). The vertical forces are due to 
elastic resistance of the deck framing or bulkheads against the motion 
of the super structure with respect to the pone The system consists, 
therefore, essentially of a beam elastically supported by another beam, 
with a shear connection to enforce equal strains at deck level. 

In this section the important simplifying assumption made is that 
the deck A-B, Figure 1, and its supports have no stiffness, and will not 
resist any relative vertical movements between hull and deckhouse. 

This simplification is not justified for any real ship ER, but 
because of its relative simplicity it is easier to study the play of 
forces; the understanding gained is of value in treating the full 
problem in the following chapters. 

In the structure shown in Figure 1, the lower hollow box beam 
represents the hull and is of length L while the upper box, the 
deckhouse, is shorter and is of length 4. Both boxes are assumed to 


be of constant cross-section. The cross-sectional area and the moment 


of inertia of the deckhouse and hull are Ay > Is and A I 


2 EA 


respectively, and the distances of the respective centers of gravity 


from each other and from the deck are a, a,a, and Q,a. 


Ih 2 


At a distance z in the free body diagram of Figure 2, the moment 


N, » and M, N, 


and direct forces in the deckhouse and hull are M, | 


respectively, with positive moments producing compression at the top of 
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the respective units. Direct forces N, and No are positive if they 
create tension. The external loads acting to the left of the section 
have a moment M. A shear force T of unknown magnitude will act on 
the underside of the deckhouse, and a similar force T will act in the 
opposite direction on the hull. Equilibrium of the portions of deck- 


house and hull in Figure 2, requires the relations: 


N OE (5a) 


De Ee (5b) 


Owing to the assumption that Navier's hypothesis is valid for the 
deckhouse and hull separately, the stresses can be determined at points 
orx 


at a distance x from the respective center of gravities. In 


1 2 
the deckhouse, - 
NM | 
1 1 
On ES (6a) 
1 Ay I, I 
in the hull, 
N M 
2 2 
A A x (6b) 
2 A, L, 2 


with tension stresses counted as positive. 

At the junction of deckhouse and hull, the longitudinal stress 
in the deck at this junction and the longitudinal stress in the deck- 
edge may differ; simple beam theory implies that they are the same. 
Shear lag in the deck plating may modify this using the following 
relation, 


0, -= ro, (3) 





soe 


where r is the ratio of longitudinal deck stress at the junction to 
the stress at deck-edge. 


Furnishing Equations (5) and (6) with xy = a Os x, =a a, and 


using Equation (3), the value of T is found to be: 


aa de tE 
see le? PAR: lL + 2 ) a 
ee i ee 1 a, I, 


T was defined as the total horizontal shear force acting between 
the left end of the deckhouse and the section at z. According to 
Equation (7), T is proportional to the moment M, and the unit horizontal 
shear (dT/dz) which will be transferred by rivets or welds from the hull 


- 


to the deckhouse, will be 


ar- oL T 


e El a 2,2 2 
u +a (a, r I, + A 1,) 


where V - SE, is the total shear force in the structure. 
After finding the value of T using Equation (7), and introducing 
this into Equations (5) and (6), the longitudinal stresses 0, and 0, 


at any point can be calculated easily, 
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CHAPTER II — GENERAL ANALYSIS OF TWO-CELL STRUCTURE 


Consider again the structure in Figure 1; differing from the 
treatment in the preceding section the assumption is made that any 
relative displacement of the deckhouse with respect to the hull will 
be resisted by the internal forces required to deflect bulkheads and 
transferse beams supporting the deckhouse. The deckhouse is considered 
as a beam elastically supported on the hull, and is further attached to 
the hull at deck level so as to enforce equal strains. External vertical 
loads and buoyancy will cause the structure to deflect, and this 
deflection can be described by the displacements Y: and Y) of the center 
lines of the deckhouse and hull respectively, (Figure 3). In order to 
exclude motions of the entire vessel as a rigid body, Yi and y, are 
defined as the relative displacements measured from a straight line 
C-D rigidly connected to the hull. As a result of this definition the 
displacement Y) of the centroid of the hull at points C and D must 


always be zero. 


Deflected 
Centerlines 








C.G. of Deckhouse 


— | nn 


C.G. of Hull FIGURE 3 
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It is assumed that the stiffness of bulkheads or deck beams, 
resisting relative vertical displacements of the deckhouse, ís constant 
for the full length of the deckhouse, the magnitude of the stiffness 
being given by a spring constant K. K is defined by Bleich as being 
the force per unit length of deckhouse required to produce a relative 
deflection equal to one unit of length. Therefore the vertical 
reaction between hull and deckhouse will be K(Yy - Yo) per unit of 
length. 

The structure analysed here is shown in Figure 4. There are two 
beams having areas Ay and A, and moments of inertia I, and Is and they 
are connected along C-D in such a way that both horizontal shear forces 
and vertical reactions can be transferred. It is assumed that, also, 


Navier's hypothesis to be valid for the hull and for the deckhouse 


separately. 


2/2 2 





loads p = = p + >», 


FIGURE 4 
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This structure will be under the action of vertical loads Py 
on the deckhouse, p, on the hull (which includes buoyancy), shear forces 


S S. and moments Mo and Mo: 


C D 
Using the "theorem of stationary potential energy", the 
differential equations for the deflections Y1 and y, Can be obtained. 


The total potential energy U consists of the internal strain energy V, 


and the potential energy pa of the external forces. The total potential 


energy U = V + U, ist, 
12 
Use [EI, y"2 + EI, y"2 + EI (a,y" +ra vie + EGM 5 
2 et 2 A” 1-1 2 1 2 
-l/2 l í 


+E ER(y'")* + E FA(y')? + E G*(y'" y™) 
1 2 a 


g 2 
>= 2pyYy 7 2p,y,] dz + [M A zes An (9) 


U will be minimum if the variation 
ôU = 0 
Using the rules of calculus of variations, the set of two simultaneous 


3 2 
equations are derived. 





See Appendix I, for the derivation of V and Ur 


> See Appendix II 
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V1 2 1V EG* V1 1V 
EE Yi + E(I, + a) I.) y] +K y] 3 + E ra, 0, I, y, 
RE (10.a) 
EG* V1 1V V1 22 1y 
enn m — * 
5 Y1 T'EST a O I, Y1 K Yi EF Y» T E(L, TAr On I) Yo 
ES Yo = Py (10.b) 


The calculus of variations method also furnish the boundary 
conditions required to determine the arbitrary constants which will 
‘appear in the general solutions of the differential equations. 


Because of symmetry, there are six boundary conditions instead of 


- 


twelve. For z = + 2/2 and z = - 2/2 
Dö 0 GLa) 
ABE y A EG e = 0 (11.b) 
EG* gr + 2EF* nn = 0 (IMIS) 
-EE* yy + E(I, + ar T y- y A N = 0 (11.d) 
EZ y + Era, 0 I, y, - EF y + E(L, + r, as I) y, =M 
(11.e) 
EEX yy - E(I, a i y + E y) - E(r 0, 0 1,) y, = 0 (11.£) 
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II.A — NON-DIMENSIONALIZATION OF THE EQUATIONS 


It is thought that employing non-dimensional equations and 


boundary conditions, the amount of algebra in the computation will 


be reduced and the results can be represented in a general form. 


Non-dimensionalization is made in the manner that is 


by the following equations. 


Ei rn 


2 
a=al 
eG 
A. = A. oe 
AA 
L 1 
P¿ =P EX, 
M=ME 2 
K=KE 


1:90 


12 


192 


i 


given 


(1222) 


(12.b) 
(12.c) 


(1274) 


(12.e) 


(12.f) 


(12.g) 


(12.h) 


The length of the deckhouse is selected as a characteristic length, 


because it is one of the most important parameters in the distribution 


of the longitudinal stresses. 


Substituting Equation (12) into Equations (10) and (11), the 


following non-dimensional equations and boundary conditions are 


derived. 


Y 
Sno 


-V1 
-a y} TaD y] 


l + a 


14 ? 


1V E o 
1 +K Y 72 Y, +b o K 
—— -y1 Ey 

SV Es U oe 


"oO 
> 


T 


ERST = 


(13.a) 


(13.b) 
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Boundary conditions for z = + 2/2 and z = - 2/2, 
Y, = 0 (14.a) 
216 Jaak + a TE = 0 (14.b) 
a an er Dn = 0 (14.0) 
“ay 6 vy. en Y nn: 7 TI Yo = 0 i (14.d) 
-a J} +by,' =a, ¥, +a, 9, =H (14.e) 
F y7 ~- aJ ta Y) “by, =0 (14.£) 


‘Expressions for non-dimensional coefficients for the equations and 


boundary conditions are given in the section entitled "Parameters 


- 


Used for Computation". 
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CHAPTER III — SOLUTION OF DIFFERENTIAL EQUATION FOR CONSTANT 
MOMENT M 


Considering the simple case that the loads Py» P, and the 


shears Sa and S_ are zero, the only loads being Ma = M, = M, 


D 
Equations (13) are then homogeneous. 

Setting the determinant of the coefficients of the differential 
equations equal to zero, the charactericits equation will be derived, 


and it will be in the following form. 


2 


4,8 6 4 _ 
r (r + a, Y + a, Y + a, Y + ay) = 0 (15) 


where, a,» ans ay» and a, are known constants. The roots of the 
characteristic equation is found for the models with different 


dimensions, and it is seen that the roots were always in the following 


manner. 
ry zr ar rS 0 (16.a) 
r =- Y (16.b) 
Lg = Y (16.c) 
r=- Ye | (16.d) 
rg = Y; (16.e) 
2 (16.£) 
ro”, +i Yo | (16.2) 
11 5 Y i Yo (16.h) 


Yio = Vy + i Yo (16.i) 
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Then, keeping in mind that the problem considered is symmetrical 
with respect to the origin of the co-ordinate z, and using only 
symmetrical functions, the general symmetrical solution will contain 


only six arbitrary constants. This general solution is, 


z 22 E = se = 
yı = Cy + Co zZz + Ca sin Y, z sinh Yı z + Cy cos Y, z cosh Y, 2 
+ Ce cosh Y, 2 = Ce cosh Y, 2 (17.a) 
- = Er = z 
Yo = Cy + C) Z + C3 sin Y, zZz sinh Ya + c, cos Y, z cosh Y Zz 
- Ha Co cosh Ya z = Uy, Ce cosh Y, 2 (17.b) 
where, 
1 = . 
C3 My Ca + H, c; (18.a) 
Pr | 
oni H3 cz + My Cy (18.b) 


and, Uys Uo: Uy, Ma are given in the following expressions, 


Cr K) -(am - bg) 
E) (-ah + bf - K) 
EE (18.c) 
(-ah + bf - K) - (am - bg) 
(am - bg) (-ah + bf - K) 


Or, 





(ah = bf + K) 


(an =- Eg) 


(- ach + a juf + K) 


(a,¿n = a,,8) 


(a 


Tee ae 


(a, gh _ auf - K) 


(-ah + bf - K) 


(am — bg) 


Or, 


(am - bg) 


(ah - bf + K) 


(- ach + a juf + K) 


(a7¿n - 27,8) 


A= 


~(a, m- a 


26 248) 


ach + a,,f + K) 
(18.d) 
248) 


(~a, ¢h + anf + K) 


-(am - bg) 


(-ah + bf - K) 
| (18.e) 
CR) 


(-ah + bf - K) 


m= a 


-(a,, 248) 


(=a,¿h + af + K) 
(SHE) 
nn 


h + a, f + K) 


24 


(18.8) 
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(- a yo + a gi +K) or (ra yê +20 ye - K) 
u 16 '4 14 4 = 4 4 
ve. A A A A, (18.h) 
(- a Ya + b mno K) (- er Y, + Ans Ya + K) 


Introduction of Equations (17) into the boundary conditions, 
leads to the equations to find the arbitrary constants given in 


r 


TABLE (1). 
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ze 


After finding the general coefficients, the stresses at any 


point along the deckhouse can be computed from the expressions for 





the moments Mj» Mos and direct forces N,» Nos 
tt 
Vt 
M, =- E L, Y, (19.b) 
Because, in this case Pj = Py "= O. 
N, =-N, = Es (a. y + 1 (20) 
1 Ae ee 


Here,it is important to notice that the solutions for y, and y, will 


be non-dimensional. 
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CHAPTER IV — SOLUTION OF DIFFERENTIAL EQUATION FOR 
- FOR EQUALLY DISTRIBUTED LOADS 


In this section the case is considered of equally distributed 
loads py and Po» acting on deckhouse and hull, respectively, while 
the moments at the end of the deckhouse are Mo = My = 0. Equilibrium 


requires external shear forces í 


2 
575,77 Cy + Py) = 


at the ends C and D. The moment in the midship section due to the 


loads pl and p2 are 


Me an T (22) 


2/2 1/2 


> z 


E 
J 


hae. 


FIGURE 5 


Equally distributed loads, 
Py FP, 





2352 


The loading being symmetrical, the general symmetrical solutions 


of Equation (13) are, 


= + i = j = = > 
Y; Cy Cy 2 Tae. SIT Y, z sinh Yi z + c, cos Yo Z cosh Yı Z 


3 
+ Cs cosh Y3 Z 4 7 cosh Ya Z 
Py + Po Be 
oe eee oe ee a 
le St ) (1 + uk 


Be. . oe a TE 1 = = 
Y, BS T Co Ze C3 sin Y, z sinh Yı z + Cy cos Y, z cosh Yı Z 


- H3 Ce cosh Y3 z -= Hy Ce cosh Ya 2 


„Pi +P? MOE O (23.b) 
24(a., +t an, + 2b) = 
14 24 (1 + HK 


1 1 . . 4 
where Cos Cy» Uy» Uy» Has and a as given in the Equations (18) for 


the case of constant moment loading, and 


5 Any, + b Vio 


But it should be mentioned that the value of "r'" shear-lag factor will 


be different in this case, so all the coefficients will not, in 
general, have the same numerical value. 

Introducing general solutions Equations (23) into the boundary 
conditions Equations (14) will lead us to the same kind of equations 


given in Table 1. But, there will be some more terms, because of the 





za 
particular integral parts of the general solutions. These additional 
terms can be determined easily, using the particular integral parts 
of the general solutions. 

After finding the general coefficients, the further computation 
to find the stresses at any point along the deckhouse follows the 
pattern for the preceding section. The only difference being that 
the additional terms appear in the equations for M, and M, due to 


distributed loads P} and p, on the deckhouse and hull, respectively. 








M, = Er g'e 25 
ley nee (25) 
lw 
tt Po i 
M, =- E L, ( > EG ) - (26) 
2W 


N, =-N RO u je 
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CHAPTER V - TOTAL SOLUTION OF DIFFERENTIAL EQUATION 


Considering the total loading for the system shown in Figure (5), 
being equally distributed loads Py and Po» acting on deckhouse and hull 
respectively, while the moments at the end of deckhouse are Mo = My = M, 
it is possible to find the total solution 

The differential equations for this system are Equations (13) 
and the boundary conditions are given by Equations (4) for +2) and 
(- =). The problem considered is symmetrical with respect to the 
origin of the co-ordinate z. The general symmetrical solutions of 
Equations (13) can be given by Equations (23.a) and (23.b) for yy and 


Yo» respectively. The coefficients in the solutions (cz, Che Hs Has 
Hz > y > Hg) are given by Equations (18) and (24). l 

But, it must be a that the values of the general 
coefficients GE Cos Cas Cys Cos Ce) will be different numerically 
than the values found for the two cases considered before. 

In the above approach to get the total solution, an additional 
assumption is made: Deck shear-lag factor "r" is considered to be 
constant over the length of the deckhouse and is given by Equation (4), 
even though there is an applied constant bending moment at the end of 
deckhouse. 

Another approach to get the total solution is to make super- 
position to the solutions found for.constant moment and equally 
distributed loading cases. 

Total solution of differential equation is found by using both 


approachs explained above, and the results are compared in the 





tr 
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CHAPTER VI — MODEL USED FOR COMPUTATION 


For this analysis a model was selected with the dimensions as 
shown in Figure (6). This model can be assumed as short deckhouse, so it 
will be possible to see more pronounced shear effects. Bulkheads are 
placed at equally spaced distances of 20 feet in the hull section. The 
thickness of the hull box girder plating is 0.5 nenes and the thickness 
of the deckhouse plating and bulkheads is 0.25 inches. The material 
constants include a Young's modulus of 30 x 10°, modulus of elasticity 
in shear of 11.5 x on and a Poisson's ratio of 0.3. 

The model is assumed to have a 15 foot draft with a corresponding 
hydrostatic upward distributed force of 17.143 tons/ft. The internal 
loading is arranged so as to provide equilibrium and a ee 
symmetrical loading. Shear and moment diagrams for the total model are 
also provided in Figure (7). 

In a model with bulkheads, the main difficulty is the determina- 
tion of the deck stiffness or spring constant (K). Since K was defined 
as the force per unit length of deckhouse required to produce a relative 
deflection equal to one unit of length, it is apparent that the value 
of K will, in reality, vary along the deckhouse length due to the 
emplacement of structural bulkheads. In order to simplify the use of 
the method, however, an average value of K must be determined. To 
achieve this end, the same approach: that was used in Reference (3) is 
followed. A symmetrical portion of the hull structure is modelled to 


include a bulkhead and attached deck and bottom plating, Figure (8). 
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SYMMETRICAL STRUCTURE FOR EVALUATION OF K 


FIGURE 8 
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This I-beam type structure is simply ee on its end and allowed 
to deflect under vertical forces applied to the hull-deckhouse con- 
nections. In the analytical procedure used for determining the forces 
needed to deflect the hull-deckhouse connections 1 inch, it must be 
kept in mind that although the deflection due to shear forces is 
negligible in most cases, in short deep metal beams the deflection 
caused by shear may become a significant portion of the total 
deflection. In this case shear contributes a major portion to the 
total deflection. 

This analytic approach to K yielded a value of 1.58 x 10° psi 
for the present model under consideration. Sample calculations are 
provided in Appendix III. | z 

In Reference (3),as a further check on K, a STRUDL program using 
'PSR' elements on the same model presented in Figure (8) was run. 
Arbitrary forces (F) were applied at the locations indicated, and the 
resulting deflection at the hull-deckhouse connection noted. The 
force was then scaled for a deflection of 1 inch. P was obtained by 
dividing F by the width of the flange (240 inches). Using Equation 
K = 2P, K was found to have value of 1.53 x 10° psi. As it is said 
in Reference (3), the disparity between the STRUDL K and the analytical 
K can be attributed to the fact that in the analytical approach the 
deflection calculations apply m the neutral axis of the beam only. 


The value of K that is used in the computations was the value 


found by employing STRUDL program. : 
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VI.A - PRESENTATION OF THE RESULTS 


Presentation of the stress distributions for Bleich's method 
and for the method developed are shown in Figures (9), (10), (11), (12), 
(13), (14), and (15) on the following pages. Figures (9) and (10) 
show the comparison of longitudinal stress distribution for constant 
moment loading at midship and at the end of deckhouse, respectively. 
The comparison of the results for equally distributed loading are shown 
in Figures (11) and (12) at midship and at the end of the deckhouse, 
respectively. Figures (13), (14), and (15) show the comparison of 
stress ee: for the total solutions. 

For the computation of the results, the model shown in Figure (6) 
is used for both methods. - 

"ACCESS II" Primer Operations in Linear Algebra for the Interdata 
Computer in Joint Computer Facility is used to find the roots of the 
characteristic equations and to solve the general coefficient matrices. 

The explanation about Bleich's method can be found in References 


Wand (3). 
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VI.B — COMPARISON OF THE METHODS 


In Figure (9), the longitudinal stress distributions at midship 
found by using two methods were plotted. In this figure, it can be 
easily seen that both results are almost the same except 400 psi 
difference in stresses at deckhouse top. The difference is 50 psi at 
main deck level. l 

Figure (10) shows the comparison of two methods at the end of 
deckhouse. In this case, the difference in stresses at deck level is 
320 psi. At deckhouse top, there is a big gap between the values of 
stresses found by using two methods. At hull bottom, the difference 
is about 30 psi. But, the value of stress found by using the theory 
developed is smaller than the value found by Bleich's eenen 

Figures (11) and (12) show the comparison of longitudinal 
stresses for equally distributed load. Looking at the results shown 
by solid line in Figure (11), shear-lag effect at main deck level can 
be seen easily. All the values for stresses are smaller than the 
results found by Bleich's method. But, the differences are not so big 
at mid ship. 

In Figure (12), it is not difficult to recognize the big 
differences between the results found by using both methods. The 
differences are approximately 1290 psi, at deckhouse top, and 690 
psi at hull bottom. At main deck level, there is about 310 psi and 
400 psi difference between the results at hullside and at the junction 


- 


of main deck and deckhouse, respectively. 
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Figure (13) shows the comparison of the total solution at 
midship. The strudl results were taken from Reference (3). The 
agreement among the results is reasonable. The stress distributions 
were transformed into moments and checked against equilibrium condition; 
that is, the values of the moments obtained corresponded to the value 
on the bending moment diagram, except for strudl results and for the 
total solution using r = 0.76, strudl results give 9% bigger than the 
moment at midship, and the total solution with r = 0.76 gives 6% smaller 
bending moment. This was expected, because shear-lag factor r = 0./6 
ie, Hot true for the whole system with applied constant bending moment 
at the end of deckhouse. 

In Figure (14), the total solutions found by using r a 0:76; 

r = 1.0 and the total solution by superposition were compared. In 
this eure, it is seen that the total solution found by using r = 1.0 
takes the average of the stresses at deck level. 

Comparison of the total solutions at the end of deckhouse were 
shown in Figure (15). The values of the stresses found by using the 
theory developed were always smaller than Bleich's results. 

As a conclusion of the comparison of the methods, inclusion of 
the shear effects into two beam theory did not change the values of 
the stresses at midship. But, at the end of deckhouse, more 


distinguishable differences were found. 
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CONCLUSIONS AND RECOMMENDATIONS 


Even though the method which includes shear effects requires 
more elaborate work, it is seen that it is possible to use it for 
computing the longitudinal stresses at any point a the deckhouse. 
There seems to be no reason to ignore these effects in design procedures 
for short and moderate deckhouses where the shear Ber could be 
pronounced. 

| After examining the comparison of the results found by using 
both methods at midship, it is possible to conclude that there is some 
indication to use Bleich's method for design purposes in the computa- 
tion of longitudinal stresses at midship. The simplicity in his compu- 
tation is a prime factor. But, for the stress solution at e end of 
the deckhouse where more distinguishable differences were found by 
employing both methods, the inclusion of shear effects into two-beam 
theory being used is necessary. The results may become more realistic. 

In the theory developed, the shear-lag effects in the deck 
represented by "r" is important only in the case of equally distributed 
load. Application of Equation (4) is a very simplified means of 
estimating "r", because it was determined by a box-girder analysis of 
the hull alone in Reference (2). Even though it may require a much 
more elaborate analysis, it is recommended to improve on the estimation 
of the shear-lag parameter. 


In this analysis, the stiffness of bulkheads or deck beams 


resisting relative vertical displacements of the deckhouse was assumed 
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constant for the full-length of the deckhouse. It is known that the 
deck stiffness K vary along the deckhouse length due to the presence 
of structural transverse bulkheads. The inclusion of this variation 
in K in the two-beam deckhouse theory needs further investigation. 
Finally, it is recommended to built a physical model similar es 
the one considered in this theoretical analysis. Then the strain and 
deflection measurement under similar loading conditions could be 


employed for verification of the results. 
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APPENDIX = I 


I.A - DERIVATION OF STRAIN ENERGY OF STRUCTURE 


Denoting by Es and E, the average longitudinal strain in the 
deckhouse and hull, respectively, the strain energy of the longitudinal 
stresses will be as follows: 


Deckhouse: 
L2 
112 


E 2 
z| (A, €, +1, Y,'0) dz 
K2 


Hull: 
2/2 
E 2 = 
> | (A, E, + L, 2 ) dz 
2/2 S 


The strains are counted positive if they represent elongation. 

In addition to the strain energy of the longitudinal stresses, 
there will be energy stored in the bulkheads or deck beams which resist 
the relative vertical displacements of deckhouse and hull; this part 
of the strain energy can be expressed by the spring constant K in the 


form 
2/2 
LO ay J dz 
1 2 
2/2 


wir 


The strain energy due to the shear deflection caused by vertical 


shear forces will be as follows: 





Deckhouse: 
Ke 2 
* 
T a ia dz 
2 
A,G 





257 





Hull: 
e 
2 u dz 
A,6 
-8/2 


where, ka and k5 are shear deflection constants and for wide flange 


box girders, can be assumed that they are equal to 1.0. 


The total strain energy V is 





2/2 
E 2 112 2 ES _ 2 
VG [A] E, tT, yp tA E +1, yg tE TY) 
-2/2 
1 2 1.2 | 
+ DEC A, 2 + EG A, Vo] dz (a) 


~ 


The stresses in the deckhouse and hull can be expressed by 
the average strains Ey and Eo» and by the second derivatives of the 


deflections ya and y 


Deckhouse: O Pen + B Ven x 


1 JE 1 


Ee FE yo. X 


Hull: O, > 


2 


Using the relationship given in Equation (3) in Chapter I 
for the longitudinal stress in the deck at the junction, and the 


longitudinal stress in the deck-edge, 
€E, =r E€ =a y + An Ya. (b) 
1 Z7 & = at) y} ra Oo. Yo 
Further, the longitudinal resultant of all stresses in the deckhouse 


NE must be equal to EA, Es and, similarly, Ny = EA, Eo: The 
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resultant of all longitudinal forces, N, + Nos must vanish. 


1 


N, + N, = ECA, E, FA, €) = 0 (c) 
By means of Equations (a) and (b), Es and E, can be expressed as 
follows: 
a E (a a + a a | 
1” TA, +A, a Gee) (d) 
aA 11 tt 
2 eA, Ee re > ©) 


Substituting expressions for €, and Ens and v, and V5 (Equations 


1 1 


1, and m, in I-C) into the Equation (a) for total strain energy V, 


~ 


2/2 
112 112 11.2 2 


tt K 
ves MY i Ey) yy 
-2/2 


t112 t112 trn ttr 


+ EX yy + Fx Y, + 6% Yi Y ] dz GE) 
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1.B - POTENTIAL ENERGY U, OF EXTERNAL FORCES 


Counting Py and P3 positive if acting downward, their potential 


energy is, 


RI? | 


-2/2 
The shear forces S and S_ and the moments M and a act 
E D C 


immediately outside points C and D, and their potential energy will 
depend on the vertical displacements Yoc and Yop? and on the rotations 
of: the end surfaces of the hull, Yoo: and Yp: 

MES into account the direction of. the shears and moments 


shown in Figure 6 in Chapter II, the potential energy will be; 


ae 


L? 2/2 
Pe ae {vee ' 
DD DO A oe lSy217 72 + My2ly y2 
and the potential energy BE of the external load is; 
2/2 
_ ei 2/2 _ 
Ihn all (pl y} + p2 y,) dz (8) 


-2/2 





SO 
I.C — EXPRESSION FOR LONGITUDINAL FORCES AND 
HORIZONTAL AND VERTICAL SHEARS 


The resultants N, and N, of the longitudinal stresses in the 


1 2 
deckhouse and hull are N, = EA. Ez > and Ny = EA, Ey» and using 


| 


Equations (d) and (e) 


EI . ; 
N, = — (a y + roa y.' 
1 a qi 2 (h) 
EL, a e 
N, = = Tr @ y] +r a, Yo ) (i) 
À. A 
2a? 
mere, L =a == > 
i A (r A, + A) 


T being the total shear force from the left end of the deckhouse to a 


point having the co-ordinate z, equilibrium requires, 
El rt tt 
A ee on) 4) 


and the unit horizontal shear will be, 


EI 


die A EE EE 
Den a e) da 
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dM 





FIGURE A 


To obtain the expression for the vertical shear Vis in the 


deckhouse, considering an element of the deckhouse of length dz, as in 


above Figure A. 


Equilibrium of moments with respect to the centroid requires, 





Substituting = from the above Equation (k) and using 


M. =- EL 


ri? 111 


on 2 
Oe a En wi en (1) 





ah 


The vertical shear force in the hull, considering Figure B below, 


similarly: 


3% UU y 


2 
Vy = - 0, 0, EI, in E(I, +r a, I.) y, (m) 


dM 


2 
A cle 





FIGURE B 


In the derivation of the Equations (1) and (m), attention must 


dM 


when the derivative of = (for i= 


1,2) is taken, the contribution from the second term in the moment 


be paid to the following point; 


equation will be zero in both cases. Because, in the following 


chapters, the cases considered for solving the equations are constant : 


moment loading (p 
d d 
mi e = 0) 
dz dz ' 


D” 0), and equally distributed loading 
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APPENDIX II - APPLICATION OF CALCULUS OF VARIATIONS 


The problem considered was in the following manner, 
2/2 ta ttt | 
U=6 F(z, Y]? Yi» Y1 » y 5 Yo: Yo: yon) O dz 


2 


ya ya TE 


To get the set of two simultaneous differential equations, 


Euler-lagrange equation is employed for y, and y., in the following 
: E 2 


manner, 
2 3 
1 
Pe u na 
dx 2 3 
dx dx 
where, - 
OF OF 
Fy = Oy? Fy = o ona 


To get the natural boundary condition equations, the following 


equations are used for Yi and Yo 


tt 


OF 
n syd z 0 


d_ 
dx 


a 20 


OF 
al ul Gym 


2 
AF d OF d OF _ 
nig gen m a) = 0 
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APPENDIX III - CALCULATIONS FOR DETERMINING K 


The following is an analytical appraoch to the determination 


of the spring constant (K) for a model with bulkheads spaced 20 feet 


Apart. 
Basic Nomenclature 
A Area of beam cross section 
K' Factor depending on shape of beam cross section 
p Distributed load at hull-deckhouse connection 
V Vertical shear due to actual forces 
V Vertical shear due to load of one pound acting at the séction 
where the deflection is to be determined i 
Yu Deflection due to internal moments 
Y Deflection due to shear 
Yr Total deflection 
X Distance along length of beam 
<> Indicates singularity functions 


The equation for the deflection due to the internal moments 
Y? is calculated through the use of singularity functions. The 
expression for the deflection due to shear e is obtained through 
the use of the method of unit loads as described in Reference (11). 


The total deflection (Yi) is then expressed as the sum of Y,, and Y 


M 
and set equal to 1 inch. The equation is solved for p (the distributed 


load acting on the hull-deckhouse connections). The expression for 


K is twice the value of p. 
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Deflection Due to Moment Only 


2 





d“y 
EI —" = 240 [p_ - p <x - 120>" - p <x - 360>7] 
x 
dx 
3 3 3 
á x? <x - 1202? _ <x - 360> 
ElYy = 240 7 er 6 + Cy x + Cy] 


B.C Y,=0 when x=0, 480" 


M 
6 
E = 30 x 10. psi 
I = 8,748,005 ae 
Y, == (2.105 x 10 ")p at x= 120" 
Deflection Due to Shear 
V | 
Y =- z —= dx (method of unit loads) 
T K' j AG 7 
y =- 1 _p(240)x 
T K“ AG 


i O E 
A O + een ren y 


nn A +m) +5 E + 9m) 


where 





I-BEAM CROSS SECTION 
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te = 055: re 0.25"; h = 3600; b= 2401 

m = 2.66; n= 0.666 

K' = 0.262 

It is customary to assume that only webs of structural shapes such as 
channels and I beams resist shearing stresses because shear stresses 


are small in flanges. 


DA 


us 
= = " 
Yr Yu =p Y 1 
eR o 
De 10° 
K = 2p K 
4 : 
K = (1.58 x 10 )psi 
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APPENDIX IV — SAMPLE CALCULATIONS 


IV.A — FOR CONSTANT MOMENT LOADING 


a = 306", A, = 150 en A, = 840 a 
I, = 534,600 in’, I, = 19,440,000 in’, 
me 31,891 900 mi 


1, = 11,917,309 a 


= A , 


M = Me = 19.2 x 107 in-lbs 


fae. 53: 108 

a, = 0-4118; a, = 0.5882, r = 1.0 (shear-lag factor) 

| Ee = =; 

De (10392) 00h, a,, = (1.232) x10, ay, = (1.136) x 10 

a = (5.685) x 10, a, = (2.335) x 10", a, = (3.129) x 10! 
4 8 6 4 2 E 

T + ay Y + a, Y + a, Y -H ag? = 0 

a, = -(1.208) x 10%, a, = (2.959) x 10%, a, = -(5.630) x 104 

ay = (1.924) x 10° 


Roots of the characteristic equation are as follows: 


4 
Br Ee 7 Ya O OS 
ee O 9.656 
B YT i Der 32200. 303.731 
ro" Yi tt Monts 3.266 +i 3.731 


11 7 Yı - i Ve 9% 2667 153.731 


= Ya teed Yo = 3.266 + i 3.731 
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f = -(583:35), g = (158.58), h= 5762.89; m= 13,700.88 


er, SS ISA MSM 
Ss, = 2236.18, S, = 2694.43 
, = (5) 0.20974, u, = (-) 0.00449, u, = (-) 0.88467 
Y, = 0.22336 
E -6 = -7 -7 

E MOS UA = COME EUR, A 
2 = 7:359 3 ome e O: Tone. a OO O 
es Nome te = 24.899 x Tore nne Toms 

_ -5 
Eg = 9-302 x 10 

-8 -8 -8 

v= 7.364 x 10 , v, = 1.065 x 19, Ya 3053/53210 
a (-)(1.305) x one Va 2.463 x DA Ve = 0.92173 tom. 
v, = 11.444 x Lon Yin ale x Ton Vy = 2.624 x 10° 

a -6 i 
EP 2755x 10 


For the B.C's equations; (z = 0,5) 


sin Yoz sinh y,z = 2.34891 
cos YZ cosh y,? = (-) 0.76987 
sin Yo? cosh y? = 2.53634 
cos Yoz sinh YZ = (-) 0.71298 
sinh Y32 ON 
cosh Y32 = 9,404 
sinh Yaz = 62.540 
z = 62.548 


cosh Y,? 
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Substituting everything into the B C's equations shown in Table 1, 


the following matrix is derived, 


1.0 0.250 -0.49611 0.1509 8.3199 -13.970||c,| = |0.0 
0.0 0.0  -3.5181 -1.049 136.010 599.812||c, 0.0 
0.0 0.0 2.549 1.297 616.40 -735.20 le 0.0 
0.0 5.249 -22,08 -94.45 -16.22 -418.6 c, 0.0 
BNO 25.51 11.25 56.73 16.12 416.3 = -3.704 x10 
0.0 0.0 37.62 14.88 9.321 769.2 6 0.0 


Results from the solution of above matrix, 


c} = 3.494 x 10°” i 
e, IAE: dor 

2 = 3,070 x 10°6 

c, == 7.899 x on 

eg = 7.021 x 10° 

eg = 2.549 x 10° 


The value of Cy is not used in the calculations, as it was pointed 


out by Bleich in Reference (1); it describes only a rigid-body motion 


of the structure not required for the purpose of this analysis. 


3 





For z = 0 (Midship) 
i 
— = - 1481.439 psi 
A 

1 
Bn 
IT Xx“ (-) 494.586 psi; x =- 126.01" (main-deck level) 

1 
del 
— x, = + 211.948 psi; x, = 54.0" (deckhouse top), 

I, i 1 

N M 
1 l 
Sos using O ===- — x 
1 As I, 1 

97 = — 1481.439 + 494.586 = - 986.85 psi (main-deck level) 
gr- 1481.439 - 211.948 = - 1693.38 psi (deckhouse top) 
D 
= = 264.542 psi 
A 

2 E 
M, 11 
i 1251.396 psi; X, = 179.989 (main-deck level) 

2 
M, n 
TT 1251.396 psi; Zn 179.989 (hull-bottom) 

2 

N M 
2 2 
Using, 0, =$ == - =— X 
2 A, I, 2 

O, = 264.542 - 1251.396 = - 986.54 psi (main-deck level) 
O, = 264.542 + 1251.396 = 1515.938 psi (hull-bottom) 


=e 





-J1l- 


For z = 0.5 (At the end of deckhouse) 


“TL 
— = ~ 982.177 psi 
A 
1 
M u 
— x, = 348.770 psi; x, = - 126.01 (main-deck level) 
I, 1 il 
M ? , 
— x, = - 149.461 psi; x, = 54" (deckhouse top) 
I 1 1 
9, = - 982.177 - 348.770 = - 1330.94 psi (main deck level) 
0, = - 982,177 + 149.461 = - 832.71 psi (deckhouse top) 
N, 
— = 175.388 psi 
A 
2 
M, 1 ES 
TX)” 1506.336 psi; X, = 179.989 (main-deck level) 
Au | 
M, n 
x, = - 1506.336 psi; xy n 179.989 (hull-bottom) 
2 


Br = 175.388 - 1506.336 


— 1330.948 psi (main deck level) 


am 175.388 + 1506.336 


1681.724 psi (hull-bottom) 





Sn 


IV.B — FOR EQUALLY DISTRIBUTED LOADING 


All the values of parameter for mathematical model will be 
the same except, 
r = 0.76 (shear-lag factor) 


A, A 


lee? _ TNG 
ie a a a 10 a 
2 
De nee in 
10 A a 


6 


1.09790 x 10”, 5 


- 1.26918 x 10 , a 2 


E 
ll 


1.09432 x 10 


214 24 


8 


a 10 a 9 


= 2,48364 x 10 `”, a y 


2.82865 x 10 


w 
il 


16 26 


Coefficients in the characteristic equation, 


ee a, = 2.98730 x 10° 


- 4.96456 x 10°, a, = 1.730691 x 10° 


= 


=] 


Roots of the characteristic equation are as follows. 


ly TY) =T¿=TY, = 0 

En Va la DS 
BE 9.624 

L ESN i y, =-~- 3.163 A i 3.621 
Yi9 = 7 ¥, + i Moen = sales t 1 2,621 
11 7 Y i ee anlo, 2723.07] 


Nv LES + 1 3.621 
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f = - 522.480, 00 2 137.393, 5 =, 700148 


m = 11,634.397 


n, = - 2.980, n, = 23.051 
r} =- 92.955, 2, = - 62.581 
S} = 1895.044, S, = 2329.223 
u, == 0.20766, u, =- 0.01375, u, =- 0.72377 
Hy = 0.25226 

_ T _ = =F 
t} =- 1.289 x 107, t,=- 8.006 x10’, t, = 3.085 x 10 

E z i E 
t, = TOEF x 10 = to = 3.697 x 10 , te = - 26.240 x 10 

26 BG a 
SNG x 107%, 25@26.307 #10, , £, = 12.570 10 
t.. = 10.272 x 107 
10 . x 
-8 -8 -8 

y, = 6.305 x 10°, v, = 1.151 x 10, v, = 25.753 x 10 

i = E > i 26 
v, = -1.855 x 10°% v, = 2.063 x 107°, v, = 0.991 x 10 

a È 26 
v, = 9.018 x 10°, vg = - 1.667 x 107%, vo = 2.367 x 10 
> 

Vio = 127041 x 10 


For the B.C's equations; (z = 0.5) 


sin Y, z sinh in Z = 2.2854 
cos Y, Zz cosh Yi z = - 0.6060 
z z = 2.4832 


sin Y, 2 cosh Ya 





he 


cosh Yo z sinh Ya, z = - 0.5578 
sinh Y, z = 9.43150 
cosh Y, z = 9.48439 
sinh y, z = 61.517 
cosh Y, z = 61.526 


Substituting everything into the B C's equations shown in 


Table 1, the following matrix is derived. 


eo 0.250 -0.4829 0.0944 6.864 — -15.55||c 


1 
0.0 0.0 -3.648 -1.268 121.0 631.2 ||e, |: 
oo 020 11.75 16.48 4943.0 -10170. Ca 
0.0 4.734 -7.452 -62.61 -38.86 -427.8 J 
0.0 24.08 7.676 62.31 35.38 A lee 
0.0 0.0 36.94 18.50 a 411.7 (lee 


JL 





Results from the solution of above matrix, 


Seen ome 
Cy = -1.000 x 10° 
Se los 707° 
pee 10015 x Ton. 
eg = -6.462 x 107 
= 3.548 x 10° 
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The value of Cy is not included in the computations. As explained 
before, it describes only a rigid-body motion of the structure, not 
required for the purposes of this analysis. 


For z = 0 (Midship) 


aT 
—— = — 957.677 psi 
A 

il 
M, 1 
—x, = - 200.573 psi; x, = - 126.01 (main deck level) 
I, a 1 
= 
e = BSS psi; x. = 54" (deckhouse top) 

I, 1 il 
using O. a = at X 

> 
I Az I 1 

a=- 957.677 + 200.573 = - 759.103 psi (main deck level) 
O = =- 957.677 - 85.592 = - 1045.62 psi (deckhouse top) 
N, 
A 

2 
M,, | 1 
7 X% = 1224.315 psi; x, = 179.989 (main deck level) 

2 
M, 11 
TXT 1226.315 psi; X E- 179.989 (hull-bottom) 

2 
Me a. 

2 A, > 2 
> = 171.371 - 1224.315 = - 1052.94 psi (main deck level) 
0, = 171.371 + 1224.315 = 1395.686 psi (hull-bottom) 





For z = 0.5 (At the end 


wi 
A 

1 
M 
= x. = = 148.120 psi; 
I 1 

1 
Hi 
I %5 63.474 psi; x] 

1 
9, = -520.778 + 148.120 
9 = — 520.778 - 63.474 
> 
DE 92.996 psi 

2 
> 
mex. = 205-474 psi: x 
I 2 

2 
HM, 
7 Xx% =- 583.474 psi; 

2 
2 
NM 
2 2 

Be x 

2 A, 9 2 
O, = 92.996 + 583.474 = 


2 


102 


of deckhouse) 


: u 
xy ST 126.01 


(main deck level) 


Sn (deckhouse top) 


- 372.658 psi (main deck level) 


- 584.252 psi (deckhouse top) 


u 
= 179.989 (main deck level) 


2 


- 490.478 psi 


676.470 psi 


x. = - 179.989 (hull-bottom) 


(main-deck level) 


(hull-bottom) 





-/7- 


IV.C — FOR TOTAL LOADING 


All the values of parameter for mathematical model will be the 
same except, 
r = 0.76 (shear-lag factor) 
7,4 


I, = 1.23670 x 10° in | 


u m m 
I = I + L, + D = 32,341,600 in 
The other coefficients and parameters will be the same as for 
equally distributed loading case. 
Substituting everything into the B.C's equations, the following 


matrix will be derived 


1.0 0.250 -0.4829 0.0944 6.864 -15.15 || c, |=] 2.322 x er 
0.0 0.0 -3.648 -1.268 121.0 631.2 |! c, | |-8.066 x 10 > 
0.0 0.0 11.75 16.48 4943.0 -10170.0 || c, | |-3.487 x 10” 
0.0 4.734 -7.452 -62.21 -38.86 -427.8 || e, | |-4.539 x 10°° 
0.0 24.08 7.676 62.31 35.38 291.5 |} c, | |-6.134 x 10% 
0.0 0.0 36.94 18.50 22.73 411.7 || cy | 2.459 x 10 | 
Results from the solution of above matrix, 

e, = 9.244 x 10°” 

Spi SDDS x w 

C} = 1.234 xX 10? 

== Ai ex Tor 

Cg == 6.169 x Tome 

8 


ce = 3.878 x 10. 





¿Bs 


For z = 0 (Midship) 


Sf 
N - 2049.499 psi 
M, 1 
— x, = - 492.638 psi; x, = - 126.01 (main deck level) 
I 1 1 
M, 
=x, = + 211.113 psi; x, = 54" (deckhouse top) ' 
I] 1 il 
a 2049.499 + 492.638 = - 1556.861 psi (main deck level) 
a 2049.499 - 211.113 = - 2260.61 psi (deckhouse top) 
N, 
— = 365.982 psi 
Ez 
M, u “ 
7%” Re x, = + 179.989 (main deck level) 
2 a 
i 
I %4 =T 3414.481 psi; Nn 179.989 (hull-bottom) 
2 
0) = 365.982 — 2414.481 = - 2048.49 psi (main deck level) 
o, = 365.982 + 2414.481 = 2780.46 psi (hull-bottom) 


For z = 0.5 (At the end of deckhouse) 


æ — 1035.029 psi 


| 
ps 
It 


u y 
T 1 463.73 psi; xX, = 126.01 (main deck level) 


| 
A 
N 


T *%]) = ~ 198.725 psi; x, = 54" (deckhouse top) 





PI 


Q 


ler 


= - 1035.029 - 463.730 = - 1498,759 psi (main deck level) 
= - 1035.029 + 198.725 = - 836.30 psi (deckhouse top) 
= 184.826 psi 
tt 

x, = 2156.883 psi; X, = 179.989 (main deck level) 

=- tt 
x, =- 2156.883 psi: Xx, = 179.989 (hull-bottom) 
= 184.826 - 2156.883 = - 1972.05 psi (main deck level) 
= 184.826 + 2156.883 = 2341.709 psi (hull-bottom) 





HO 


li. 


2, 


RS 


14. 
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